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Abstract: Mid-western China is one of the most sensitive and fragile areas on the Earth. 
Evapotranspiration (ET) is a key part of hydrological cycle in these areas and is affected by both global 
climate change and human activities. The dynamic changes in ET and potential evapotranspiration (PET), 
which can reflect water consumption and demand, are still unclear, and there is a lack of predictive 
capacity on drought severity. In this study, we used global MODIS (moderate-resolution imaging 
spectroradiometer) terrestrial ET (MOD16) products, Morlet wavelet analysis, and simple linear regression 
to investigate the spatiotemporal variations of ET, PET, reference ET (ET»), and aridity index (AI) in 
mid-western pastoral regions of China (including Gansu Province, Qinghai Province, Ningxia Hui 
Autonomous Region, and part of Inner Mongolia Autonomous Region) from 2001 to 2016. The results 
showed that the overall ET gradually increased from east to southwest in the study area. Actual ET 
showed an increasing trend, whereas PET tended to decrease from 2001 to 2016. The change in ET was 
affected by vegetation types. During the study period, the average annual ET» and AI tended to decrease. 
At the monthly scale within a year, AI value decreased from January to July and then increased. The 
interannual variations of ET) and AI showed periodicity with a main period of 14 a, and two other 
periodicities of 11 and 5 a. This study showed that in recent years, drought in these pastoral regions of 
mid-western China has been alleviated. Therefore, it is foreseeable that the demand for irrigation water for 
agricultural production in these regions will decrease. 
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1 Introduction 


Mid-western China, an arid and semi-arid region, is ecologically important and plays 
an important role in supporting the livestock husbandry in China. It is one of the most sensitive 
and fragile areas on the Earth, yet readily influenced by global climate change and human 
activities. The increase in greenhouse gas concentrations due to human activities has led to 
surface warming over almost the entire globe (Stocker et al., 2013). Concomitantly, the 
climate in mid-western China has also been shifting from warm-dry to warm-wet (Shi et al., 
2006). Even minor variations in precipitation and temperature can induce significant changes in 
hydrological processes such as 
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evapotranspiration (ET) (Gan, 2000; Ma et al., 2004). Additionally, the land surface water cycle 
and ET are affected by global climate change (Jung et al., 2010; Meng, 2021). Therefore, further 
documenting the changes in ET and linking them with variations in vegetation cover will 
contribute to a more comprehensive understanding of agricultural irrigation and ecological water 
demand, thus promoting the development of agriculture and livestock husbandry. 

ET includes the evaporation of water from the Earth's surface and the transpiration of 
vegetation (Chahine, 1992). It returns about 60.0% of the global terrestrial precipitation to the 
atmosphere (Wang and Dickinson, 2012) and consumes about 60.0% of the net surface radiation 
energy (Trenberth et al., 2009). ET links the terrestrial water cycle and surface energy 
transmission, which determines the intensity of the interactions between the Earth and atmosphere 
(Jung et al., 2010). In arid areas, global climate change may increase (or decrease) potential 
evapotranspiration (PET), thereby increasing (or decreasing) the aridity (Huo et al., 2013). 
Human activities, such as domestic animal husbandry, also play a role in influencing water 
cycling (Goyal, 2004). In recent years, the potential impact of climate change on reference 
evapotranspiration (ETo) has been extensively studied (Wang et al., 2019; Liu et al., 2020; Jerin et 
al., 2021). Long-term climate sequence analysis showed that abrupt climate change occurred in 
northwestern China around 1975 (Chen et al., 1991). Time series analysis revealed that ETo 
decreased in all seasons from 1954 to 1993 in China, especially in northwestern and southeastern 
areas (Wang et al., 2007; Zhang et al., 2007; Song et al., 2010; Zhang et al., 2010; Thomas, 2015). 

Some researchers have analyzed the factors influencing ET changes, revealing that the 
sensitivity of ETo changed according to location and associated environmental variables (Liu et 
al., 2010). For example, Liu et al. (2010) reported that temperature increase was the main cause of 
ETo rise in the Yellow River Basin of China. Eslamian et al. (2011) reported that temperature and 
relative humidity were two most important variables affecting the monthly trend of ETo in Iran. 
Espadafor et al. (2011) pointed out that the increases of temperature and solar radiation and the 
decrease of relative humidity were the main factors leading to the increase of ET in southern 
Spain. Dinpashoh et al. (2011) proved that wind speed was the main cause of ET» rise in Iran. Li 
et al. (2013) found that the influence of meteorological factors on ETo was different in the Heihe 
River Basin of China at the spatiotemporal scales. Mosaedi et al. (2016) found that sunshine 
hours, wind speed, and temperature had great effects on ETo in Iran. Since the late 1980s, some 
national construction projects and the adjustment of agricultural planting approaches have been 
implemented in mid-western China, which has affected the environment, including surface 
vegetation cover. Further exploring ET change and its links with vegetation cover change would 
be of great help to optimize water management for agriculture production and livestock 
husbandry in mid-western China. 

Ground-based observation can provide data to show time series ET information, but it is 
difficult to exhibit the spatial distribution characteristics of ET for a broad area when there are not 
enough observation sites, for instance, in arid mid-western China. By contrast, remote sensing 
technology can help to obtain the spatial variability of ET at a regional or global scale. Based on 
the Penman-Monteith equation and moderate-resolution imaging spectroradiometer (MODIS), the 
global MODIS terrestrial ET (MOD16) data have been released with an accuracy of 86.0% (Mu et 
al., 2007, 2011). The MOD16 products have been validated and applied in many areas with 
various climate conditions (Jia et al., 2012; Yang et al., 2014; Chen et al., 2015; Sun et al., 2021). 
Jia et al. (2012) found that the validation of remotely sensed-ET (RS-ET) at the basin scale 
showed a good consistency between the 1-km annual RS-ET and the validated data, such as the 
water balance ET, MODIS ET products, precipitation, and land-use types. Kim et al. (2012) 
reported that the MOD16 products performed best at five forest sites and were mismatched with 
observed ET data at two grassland sites in Asia. Feng et al. (2015) pointed out that the MOD16 
algorithm reduced the accuracy of ET estimations for grassland, savanna, and shrubland sites over 
semi-arid ecosystems. In the Tibetan Plateau of China, Chang et al. (2018) also obtained good 
consistency between the MOD 16 products and the measured values. However, the performance of 
the MOD 16 products varies under different climate and underlying surface conditions. Therefore, 
further investigation on spatiotemporal characteristics of ET using the MOD16 products in 
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mid-western China would aid understanding of the spatiotemporal distributions and variations of 
regional ET under global climate change and human activities. 

In this study, the MOD 16 products were used to analyze the spatiotemporal changes of ET and 
PET, and ground-based data were used to explore the changes in ETo and aridity index (AI) in 
typical mid-western pastoral regions of China, including Gansu Province, Qinghai Province, 
Ningxia Hui Autonomous Region, and part of Inner Mongolia Autonomous Region. The 
objectives of this study were to: (1) analyze the spatiotemporal variations of ET and PET; (2) 
reveal the temporal variations of ET with different vegetation types; and (3) explore the temporal 
variations of ETo and AI. It is expected that this study will provide a reference for understanding 
agricultural irrigation and ecological water demand in mid-western China. 


2 Materials and methods 


2.1 Study area 


This study was conducted in typical mid-western pastoral regions of China, including Gansu 
Province, Qinghai Province, Ningxia Hui Autonomous Region, and part of Inner Mongolia 
Autonomous Region. The longitude range is 89°30’-113°14’E and the latitude range is 33°02’— 
43°28'N. The elevation is 562—6820 m and the terrain gradually decreases from southwest to 
northeast. These pastoral regions are characterized by an arid and semi-arid climate with high 
radiation and heat resources, little precipitation, strong evaporation, and large temperature 
difference between day and night (Meng et al., 2013; Liu et al., 2021). 


2.2 MOD16 products 


The MODIS ET dataset includes surface ET, latent heat flux, PET, and potential latent heat flux. 
The dataset has a spatial resolution of 1 kmx1 km and time resolutions of 8 d, 1 month, and 1 a. 
The algorithm is based on the Penman-Monteith equation with some improvements (Mu et al., 
2007). These data can better reflect the non-uniformity of the underlying surfaces of desert and 
oasis, and are suitable for the study of surface ET in arid areas. The annual synthetic product 
(MOD 16) with actual surface ET data was selected, and the time series was from January 2001 to 
December 2016. The original MODIS products were stored in hierarchical data format (HDF) 
with sinusoidal projection. HDF files of the MOD16 products were converted into GeoTIFF files 
using ArcGIS 10.3 and ENVI+IDL (The Environment for Visualizing Images+Interactive Data 
Language) software. The projection conversion, raster mosaic, resampling, and other operations 
were performed to organize the final ET data. 


2.3 Vegetation and meteorological data 


The vegetation dataset was provided by the Data Center for Resources and Environmental 
Sciences, Chinese Academy of Sciences (RESDC) (http://www.resdc.cn). This dataset includes 
vegetation types, horizontal and vertical distributions, and the relationship of vegetation with 
climatic factors and surface environmental factors. Distribution of vegetation types in the study 
area is shown in Figure 1. 

The meteorological data (including daily evaporation, precipitation, average temperature, 
average pressure, relative humidity, and average wind speed) were obtained from 94 
meteorological stations (Fig. 1) in the study area from 2001 to 2016. For comparative analysis 
with ET data, the daily values were processed into 192 monthly and 16 annual values. 


2.4 Reference evapotranspiration (ETo) and aridity index (AI) 


ETo reflects an important component of the atmosphere-hydrosphere water cycle, which drives 
the dry-wet conditions on the Earth's surface. The Penman-Monteith equation was recommended 
by the Food and Agriculture Organization (FAO) for the calculation of ETo (Allen et al., 2007); it 
was adopted in this study. The formula below was used to calculate ETo: 
_ 0.408A(R, — G) + 900y / (T +273)U (e, —e,) 
A+y(1+0.34U,) j 


ETy (1) 
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Fig. 1 Distribution of vegetation types and location of meteorological stations in the study area 


where ETo is the reference evapotranspiration (mm); Rn is the net radiation at the crop surface 
(MJ/(m?-d)); G is the soil heat flux density (MJ/(m?.d)); T is the mean daily temperature at 2 m 
height (°C); U2 is the wind speed at 2 m height (m/s); es is the saturation vapor pressure (kPa); ea 
is the actual vapor pressure (kPa); A is the slope of the vapor pressure curve (kPa/°C); and y is the 
psychometric constant (kPa/°C). 

To further analyze dry-wet conditions in terms of agriculture demand, we also calculated AI. In 
arid areas, the agricultural aridity is induced mainly by ET exceeding precipitation. ETp means the 
reference crop ET determined by climate factors, which can be easily observed (Huo et al., 2013). 
AI, defined by Thornthwaite (1948), represents the degree of drought in arid or semi-arid areas 
and is calculated as: 


a (2) 


where AI is the aridity index and P is the precipitation (mm). AI can also show the water required 
to satisfy ETo. If there is no precipitation, AI is equal to 1, and the mean aridity is the highest. By 
contrast, if precipitation is equal to or higher than ETo, AI is zero or negative. In the study area, as 
ETo is greater than precipitation, AI is between 0 and 1. 


2.5 Interannual variation of evapotranspiration (ET) 


The one-way linear regression method (Feng et al., 2015) was used to analyze the temporal 
variations of annual average ET in the study area with the formula: 


n n n 
mix, LIL, 
jal j=l j=l 
2 3 (3) 
z Pa z . 
IT -| Ds 
j=l j=l 
where Osiope is the slope of the trend line; n is the accumulated years of monitoring (a); and X; is 


the annual ET of the j year (mm). If Osiope>0, ET increases during the monitored n years, and vice 
versa. 


Olope = 


2.6 Data analysis 


Morlet wavelet analysis has been widely used in many fields, showing the periodic distribution, 
amplitude, and bit equality information of time series in time domain, and has the ability to 
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diagnose mutation points in mathematical sense (Ma et al., 2013; Wu et al., 2019; Nourani et al., 
2021). Its basic principle is to use a cluster of wavelet function system to represent or 
approximate a signal or function (Burrus, 2005). 

Wavelet function (1) € L7(R) is as follows: 


[wat =o, (4) 
where ¢ is time variable (a) and Y(t) is a basic wavelet function and can form a cluster of functions: 
1 
—,(t-b 
vaak) (5) 


where Y4, (f) is the wavelet transform function; a is the expansion scale; and b is the translation 
parameter (a, b E R, a#0). 
For f(t) € L’(R), the continuous wavelet transform is as follows: 


1 — amy 
valh ror) dt, (6) 


where f(t) is a signal or square integrable function; z(=) is a complex conjugate function of 


a 
t—b 
(=) | 
a 
Wavelet variance can be used to determine the relative intensity of different scale disturbances 


in the signal and the main time scale, i.e., the main period (Burrus, 2005). The formula is as 
follows: 


io 2 
Var(a) = D w, (a, b)| db, (7) 


where Wa, b) is wavelet coefficient. 

In this study, the wavelet toolkit of MATLAB software was used for wavelet transform. In 
addition, data processing, analyzing, and drawing were carried out in ArcGIS 10.3, ENVI 5.6, 
MATLAB 2016, and EXCEL 2010 software. 


3 Results 


3.1 Validations of the MOD16 data 


The accuracy of the MOD16 data was verified using the correlation coefficient between PET 
calculated from the MOD16 products and ET observed from the meteorological stations. As shown 
in Figure 2, at the temporal scale, the correlation coefficient between calculated PET and 
observed ET was 0.9008 (Fig. 2a) and at the spatial scale, the correlation coefficient between them 
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Fig. 2 Relationships between potential evapotranspiration (PET) calculated from the MOD16 products and 
evapotranspiration (ET) observed from meteorological stations along a chronosequence (a) and among sites (b) in 
mid-western China. The data from 192 time points and 94 meteorological stations were used to analyze the 
relationship between PET calculated from the MOD16 products and ET observed from meteorological stations 
from 2001 to 2016. 
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was 0.5993 (Fig. 2b). This showed that PET calculated from the MOD16 products was consistent 
with ET observed from the meteorological stations, and this dataset from the MOD16 could be 
used to analyze the spatiotemporal distribution of land surface ET in mid-western China. 


3.2 Spatiotemporal variations of ET and potential evapotranspiration (PET) 


The overall ET in the study area gradually increased from east to southwest (Fig. 3). During the 
studied 16 a, ET varied from region to region across the study area, but PET was stable. 


PET (mm) 
mm <500 
= 500-700 
m= 700-900 
== 900-1100 
== 1100-1300 
== 1300-1400 
= 1400-1500 
= 1500-1700 
0 200 400 km m= 600-700 = 1700-1900 
—1— ma >700 m= >1900 

[= No data © No data 


(G) 


Slope of PET (mm/a) 
mm <-20 


Slope of ET (mm/a) 
mm <-20 


mm -20--10 =m -20--10 
= -10--3 mm -10--3 
= -3-3 mm -3-3 

= 3-10 m™ 3-10 

m= 10-20 m= 10-20 
me >20 ma >20 

C No data CI No data 


Fig. 3 Spatial variations of ET (a) and PET (b) as well as spatial variations of slopes of ET (c) and PET (d) in 
mid-western China. Note that due to the restriction of the MOD16 products, ET in areas without vegetation such 
as deserts was not calculated and these areas appeared to be blank (white). 


From 2001 to 2016, the average annual ET in the study area increased and the difference in 
magnitude among pastoral regions was obvious, with the order of Qinghai>Gansu>Ningxia> 
midwestern Inner Mongolia (Fig. 4). The increasing rates of ET were 0.35 mm/a in Qinghai, 8.39 
mm/a in Ningxia, 5.90 mm/a in Gansu, and 2.51 mm/a in midwestern Inner Mongolia. The average 
annual PET, with the range of 980.00—1500.00 mm, showed a decreasing trend during the study 
period. The order of PET among the four regions was midwestern Inner 
Mongolia>Ningxia>Gansu>Qinghai. The decreasing rates of PET were 7.00 mm/a in midwestern 
Inner Mongolia, 11.10 mm/a in Ningxia, 12.00 mm/a in Gansu, and 5.42 mm/a in Qinghai. The 
difference between PET and ET can explain the water shortage on the land surface, that is, the 
drought conditions. The differences between PET and ET were large in all four regions, indicating 
that the overall water condition was dry and water deficit existed. The degree of drought among the 
four regions was in the order of midwestern Inner Mongolia>Ningxia>Gansu>Qinghai. The trends 
of PET were decreasing, showing that drought conditions in different regions were mitigated. 


3.3 Average monthly ET and PET changes 


The trend of average monthly ET and PET changes in various pastoral regions was dominated by 
the higher ET in summer and the lower ET in other seasons, and both ET and PET tended to 
increase first and then decrease at the monthly scale (Fig. 5). The maximum monthly ET in 
Gansu, midwestern Inner Mongolia, Ningxia, and Qinghai was 61.85, 27.50, 36.20, and 61.10 
mm, respectively. The minimum monthly ET in Gansu, midwestern Inner Mongolia, Ningxia, and 
Qinghai was 24.63, 10.79, 16.17, and 26.32 mm, respectively, which respectively occurred in 
December, May, April, and December. The minimum monthly PET in the study area was in 
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Fig. 4 Average annual variations of ET and PET in mid-western China (a) as well as in the pastoral regions of 
Gansu (b), midwestern Inner Mongolia (c), Ningxia (d), and Qinghai (e) during 2001-2016 


December, with the values of 51.31, 34.35, 45.03, and 43.60 mm in Gansu, midwestern Inner 
Mongolia, Ningxia, and Qinghai, respectively. The maximum monthly PET was 209.03 mm in 
Gansu, 229.10 mm in midwestern Inner Mongolia, 219.27 mm in Ningxia, and 190.85 mm in 
Qinghai. The maximum difference between the average monthly PET and the average monthly 
ET in Gansu and midwestern Inner Mongolia was 153.90 and 207.80 mm in May, respectively, 
whereas it was 190.53 mm in June for Ningxia and 129.75 mm in July for Qinghai. 


3.4 Interannual variations and monthly changes of ET in different vegetation types 


From 2001 to 2016, the interannual variation of ET for each vegetation type showed a steadily 
increasing trend (Fig. 6). The average annual ET over the study period in different vegetation 
types was as followed: forest (386.85 mm), shrub (337.68 mm), alpine vegetation (303.77 mm), 
natural grassland (296.50 mm), cultivated vegetation (286.39 mm), swamp (279.19 mm), and 
desert (187.72 mm). 

The change of monthly ET was significantly different in different vegetation types (Fig. 6). The 
trend can be roughly divided into three categories. First, there was a single peak mainly for forest; 
specifically, monthly ET increased from January and decreased after reaching the maximum in 
July at the monthly scale. Second, monthly ET fluctuated in natural grassland, desert, cultivated 
vegetation, and swamp; specifically, monthly ET firstly decreased from January to April and then 
increased to peak in July and/or August. After reaching the peak value, monthly ET decreased and 
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later increased again. Third, monthly ET exhibited a bimodal type for alpine vegetation and 
shrubs. For this category, monthly ET first slowly increased till March, slightly decreased in 
April, then rapidly increased to peak in July, and finally decreased. 
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Fig. 5 Average monthly ET and PET changes in mid-western China (a) as well as in the pastoral regions of 
Gansu (b), midwestern Inner Mongolia (c), Ningxia (d), and Qinghai (e) during 2001-2016 
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Interannual variations (a) and monthly changes (b) of ET in different vegetation types in mid-western 
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3.5 Interannual variations and monthly changes of ETo, AI, and precipitation 


Linear regression analysis showed that the average annual ETo and AI in the study area decreased 
from 2001 to 2016 (Fig. 7). Specifically, the average annual ETo demonstrated a significant 
decline trend at a rate of 12.90 mm/a, and the average annual AI decreased by 0.03 mm/decade. 
By contrast, the average annual precipitation increased from 2001 to 2016, at a rate of 1.98 mm/a. 
It should be noted that the changes of average annual AI and the average annual precipitation 
were not statistically significant. 

The average monthly ETo, precipitation, and AI had significant seasonal variations (Fig. 7b). 
Overall, both ETo and precipitation increased first and then decreased at the monthly scale, while AI 
exhibited the opposite trend and was greater than 0.60. This indicated that the study area belonged 
to an arid area and the climate was dry throughout the year. In July, ETo reached the maximum with 
a value of 183.94 mm, and precipitation also had the maximum value in July with a value of 63.34 
mm. Both the temporal trends of ETo and precipitation caused a minimum AI of 0.60 in July. During 
the periods of January—April and November—December, AI was higher than 0.90. 
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Fig. 7 Interannual variations (a) and monthly changes (b) of average reference evapotranspiration (ETo), 
precipitation, and arid index (AI) in mid-western China 


3.6 Periodic features of ET» and AI 


Morlet wavelet analysis revealed that ETo and AI had obvious localization characteristics at the 
interannual time scale, with nested signal strength and phase, and uneven frequency and time 
domain distribution (Fig. 8). The variation of annual ETo had a long periodicity of 10-16 a, and 
there was also a short periodicity of 4—6 a. Therefore, annual ETo had three periodicities: 5, 11, 
and 14 a, indicating that the change of annual ETo had a major periodicity of 14 a, a secondary 
periodicity of 11 a, and a minor periodicity of 5 a. The major and secondary periodicities of 
annual AI were consistent with those of annual ETo, but annual AI still had shorter periodicities of 
about 3 and 6 a. 


4 Discussion 


4.1 Applicability of the MOD16 products in arid and semi-arid areas 


In recent years, the MOD16 products are widely used in hydrological research in arid and 
semi-arid areas where ground-based data are scarce (Khan et al., 2018; Li et al., 2021; Wang et 
al., 2021). For example, Glenn et al. (2011) used the MOD16 products to capture the general 
actual ET as well as its features and trends at a continental scale in Australia. Mu et al. (2011) 
compared the MOD16 derived ET data with measured ET data from 46 AmeriFlux eddy 
covariance flux towers located across seven biomes in North America and found that the average 
error was 24.1%, within the accuracy range of multiple remote sensing methods (10.0%—30.0%) 
(Kalma et al., 2008). Moreira et al. (2018) evaluated the performance of the MOD16 products in 
nine eddy covariance monitoring sites from the Large-Scale Biosphere-Atmosphere Experiment 
in the Amazon and revealed that the products showed similar ET values with the calculated 
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Fig. 8 Morlet wavelet analysis of periodic features of average ETo (a, b) and AI (c, d) in mid-western China 


statistics. Although the accuracy of the MOD16 products still needs to be further verified in 
semi-arid areas, it is helpful to identify the spatiotemporal trends of ET change (Jovanovic et al., 
2015). Concerns have also arisen on the short duration of the MOD16 dataset (only from the last 
20 a). Therefore, it is better to use additional datasets from the MOD16 products over a longer 
duration to demonstrate ET change associated with climate change. 


4.2 Effects of vegetation type and climate change on ET 


In this study, actual ET in mid-western China increased year to year, although the trend was not 
significant. This is because the amount of ET is closely related to both vegetation type and 
climate change. Vegetation connects the soil to the atmosphere and regulates water balance, so 
higher vegetation cover can facilitate water cycling. In recent years, global warming and 
significant increase in vegetation cover in Northwest China (Mou et al., 2018) have led to an 
increase in ET. Mo et al. (2014) showed that NDVI (Normalized Difference Vegetation Index) 
was positively correlated with ET in China. We also found PET exhibited a slightly decreasing 
trend with years in mid-western China, which is consistent with the finding in Northwest China 
revealed by Yang et al. (2015). The declining trend of PET suggested that drought stress was 
ameliorated (Yao et al., 2015). Precipitation in Northwest China was gradually increasing in 
recent years (Chen et al., 2014), and there was a significant negative relationship between 
precipitation and PET (Pan et al., 2021). Additionally, the variation of PET was possibly affected 
by other atmospheric factors, such as solar radiation, wind speed, and temperature (Baguis et al., 
2010; Thomas, 2015; Feng et al., 2019). 

In this study, ET in mid-western China changed in an obvious season-specific way. There was 
little precipitation from January to April in this area and surface water could not be replenished, 
which results in a decreasing tendency of ET. In April, vegetative plants began to grow, leading 
to an increase of ET. In the northwest pastoral areas, the land use/cover types were complex, 
including desert, wind- and sand-controlled engineering areas, and the forests and grasslands 
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were converted from farmlands. These areas tended to have a low vegetation cover and the 
existing vegetation was mainly cultivated vegetation that can absorb shallow soil moisture. 
Therefore, the increase of ET was not significant. In July and August, precipitation was 
concentrated, thus ET reached the highest value. Moreover, solar radiation was a major factor 
affecting ET (He et al., 2015). Summer days were long, receiving the most solar radiation and 
promoting the growth of natural and cultivated vegetation. Surface water circulation was 
accelerated from June to August, and ET continued to increase. In this study, the seasonal 
response of ET varied in different pastoral regions. The lowest ET in mid-western Inner 
Mongolia and Ningxia occurred in May and April, respectively, and ET in Gansu was also at a 
minimum value in April. This is probably because of the snow melting in spring and the 
gradually increased temperature. Appropriate temperature and moisture were beneficial for 
grass or crop re-greening, and ET in this period was lower than that in other times of the year 
(Tong et al., 2016). 

The ET of all individual vegetation types changed similarly to the overall ET in mid-western 
China. However, vegetation types varied in different pastoral regions of the study area and land 
use/cover types were diversified, resulting in at least three types of monthly ET change trends. 
We found there was a decrease in ET around April for some vegetation types. In spring, natural 
vegetation began to regrow, but grazing and human activities (such as plowing) would lead to a 
reduction of vegetation cover. Considering low precipitation in this season, the decrease in ET 
can be predicted. After spring, temperature and precipitation gradually increased, providing 
adequate heat and water for the rapid growth of vegetation, which leads to a rapid rise of ET to 
peak around August. 


4.3 Periodic changes of AI and ETo 


The variation of ETo was influenced by radiometric and aerodynamic variables, which mainly 
reflect the changes in topography and climate (McVicar et al., 2007; Huo et al., 2013; Collins et 
al., 2021). In this study, ETo fluctuated and decreased during the studied 16-a duration. This 
was consistent with the studies conducted in India (Chattopadhyay and Hulme, 1997), and the 
Yangtze River Basin (Xu et al., 2006) and the Yellow River Basin in China (Zhang et al., 2015). 
The decreasing trend of ETo was mainly attributed to the decreases in net total radiation and 
wind speed (Zhang et al., 2015). However, Tabari (2010) showed that average annual ETo 
increased by 2.30—11.30 mm in an arid western area of Iran, which may be related to the higher 
temperature. Nouri and Bannayan (2019) suggested that the variation trend of ET) was 
regulated by mean temperature and wind speed in humid environment. 

In this study, we also found that AI decreased in mid-western China, suggesting the whole 
area became wetter. This was consistent with previous results (Wang et al., 2007; Liu et al., 
2013). Vegetation cover has increased in this study area during the period from 2001 to 2015 
(Mot et al., 2018) and there was a humidification trend in Northwest China based on the study 
of a dry-wet homogenization index (Wang et al., 2007). Analysis of the monthly changes of ETo 
and AI showed that there was a unimodal change for ET» and an opposite change for AI within 
a year in this study. Past studies revealed that temperature played a decisive role in regulating 
water cycle in most arid areas, such as Iran (Bazrafshan, 2017). AI was significantly and 
negatively correlated with precipitation, which contributed 91.7% of the decrease in AI in 
Northwest China from 1960 to 2010 (Liu et al., 2013). In this study, temperature and 
precipitation in mid-western China were higher from June to August, leading to the greatest 
ETo; during this period, AI was the lowest, indicating the wettest time of the year. Zhang et al. 
(2015) found that the periodic periods of ETo and AI were 4—6, 8-12, and 18-22 a in the Yellow 
River Basin, China. The periodicity of ETo was similar to that of AI to some extent, and the 
major periodicities of the two variables were comparable. Also, the periodic changes of AI and 
ETo were affected by macroclimate patterns, such as El Niño, and atmospheric circulation 
(Zhang et al., 2015). Atmospheric circulation had a periodicity of 2—4 a, and it took 2—7 a for 
each El Niño to occur (Fan, 2018). However, the relationship of periodic changes of AI and ETo 
with global warming still needs to be explored with long-time monitoring. 
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5 Conclusions 


The overall ET gradually increased from east to southwest in mid-western China. From 2001 to 
2016, actual ET showed an increasing trend, whereas PET decreased. The average annual ET 
variation in different vegetation types can be divided into unimodal, fluctuation, and bimodal 
patterns. The average annual ETo and AI decreased during the studied 16-a duration. The 
interannual variations of ETo and AI showed obvious periodic changes with a major periodicity of 
14 a, a secondary periodicity of 11 a, and a minor periodicity of about 5 a. This study showed that 
drought has been alleviated in mid-western pastoral regions of China in recent years. Therefore, a 
decrease in irrigation water demand may be foreseen for agricultural production. This research 
will help clarify irrigational and ecological water requirements and guide irrigation management 
for the development of agriculture and animal husbandry in arid areas. 
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